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Abstract. We consider a system of PDEs of Monge-Kantorovich type that, in the 
isotropic case, describes the stationary configurations of two-layers models in granular 
matter theory with a general source and a general boundary data. We propose a new 
weak formulation which is consistent with the physical model and permits us to prove 
existence and uniqueness results. 



1. Introduction 

The model system usually considered for the description of the stationary configura- 
tions of sandpiles on a container is the Monge-Kantorovich type system of PDEs 

- div(v Du) = f in 0, 

Du\ < 1, v > in O, 

(1) ^(l-|D«|)u = in SI, 

u < (f) on dSl, 

u = (f> on Tf 

(see, e.g., HI HJJ IS])- The data of the problem are the flat surface of the container 
SI C M 2 , the profile of the rim (f>, and the density of the source / > 0, whereas the set 
Tf is a subset of dSl, defined in terms of the other data, that will be specified below. 

The dynamical behaviour of the granular matter is pictured by the pair (u, v), where 
u is the profile of the standing layer, whose slope has not to exceed a critical value 
(\Du\ < 1) in order to prevent avalanches, while v > is the thickness of the rolling 
layer. The condition (1 — \Du\)v = corresponds to require that the matter runs down 
only in the region where the slope of the heaps is maximal. 

The set V f (which depends on the source /, the geometry of SI, and on the boundary 
datum cj)) is the part of the border where every admissible profile u touches the rim, in 
such a way the exceeding sand can fall down (see Definitions in Section [3]). We underline 
that the set r f is not an additional datum of the problem, but it is constructed in terms 
of the other data (see ((9j) for its precise definition). 

The main contribution of our results to the theory concerns the uniqueness of the 
u-component for general boundary value problems, based on a new weak formulation of 
the continuity equation — div (v Du) = f in Si. 
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The case of the open table problem, corresponding to u = (ft = on d£l, is already 
completely understood (see e.g. |H O H [T7] and the references therein). Namely, if 
da denotes the distance function from the boundary of Q, it is possible to construct a 
function Vf > 0, Vf € L 1 ^) such that the pair (dn,Vf) is a solution to ([!]) (we underline 
that the continuity equation is understood in the sense of distributions). Moreover 
it turns out that Vf is the unique admissible f-component, and every profile u must 
coincides with where the transport is active. 

These results validate the model for the open table problem, since they depicted the 
sole physically acceptable situation: the mass transport density v has to be uniquely 
determined by the data of the problem, while the profile u could be different from the 
maximal one only where the mass transportation does not act. 

Moreover the profile is unique (and maximal) if and only if the source / pours sand 
along the ridge of the maximal profile (i.e. on the closure of the set where do, is not 
differentiable) . 

As far as we know, only the following two particular cases of non-homogeneous bound- 
ary conditions were considered in literature. 

In [13], mostly devoted to a numerical point of view, the problem of the open table 
with walls (corresponding to u = (ft = on a regular portion T of dQ, and (ft = +oo in 
dQ \ r) is considered. In order to take into account the fact that the sand can flow out 
from the table only through T, the weak formulation of the continuity equation proposed 
in [T3] is the following: 

/ v(Du, Dip)dx = [ fijtdx, V^C c °°(l 2 \r). 
Jn Jn 

Under suitable regularity assumptions on the geometry of the sandpile, it is proved that 

there exists a function vj > 0, vt £ L 1 ($7) such that the pair (dr,Vf) is a solution to 

([!]), where dr is the distance function from T. 

A different approach to non-homogeneous boundary conditions was recently proposed 

in [12]. In that paper we considered only admissible boundary data, that is continuous 

functions (ft on d£l that coincide on the boundary with the related Lax-Hopf function 

U(f,. In the model, this corresponds to treating the so called tray table problem, where 

the boundary datum (ft gives the height of the rim. The requirements are that the border 

of the rim is always reached (u = (ft on d£l), and that the continuity equation is satisfied 

in the sense of distributions. The existence of a solution is obtained, in analogy with 

the open table problem, by exhibiting an explicit function Vf > 0, vj 6 L 1 ^) such 

that the pair (u^,Vf) is a solution to ([!]). Moreover a necessary and sufficient condition 

for the uniqueness of the u-component can be obtained, with minor changes, as in the 

homogeneous case. 

The main novelty in the analysis of the non-homogeneous case concerns the lack of 
uniqueness of the ^-component. Namely, the boundary datum (ft modifies the geometry 
of the directions along which the sand falls down. In particular it may happen that 
a family of transport rays passing across f2 covers a set of positive measure, so that 
it is possible to transport any additional mass along these rays, keeping the total flux 
unchanged. 

In the present paper we shall deal with general boundary data, thus allowing the 
presence of walls as well as of exit points at different heights. The main goal will be to 
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modify the weak formulation of the continuity equation in order to gain the uniqueness 
of the u-component, without loosing information concerning the n-component, then 
validating the model in a very general case. 

Moreover we shall not require that the profiles have to reach the height of the rim 
at every point of dQ where they a-priori could agree, compatibly with their gradient 
constraint (i.e. at every point where the maximal profile ua, agrees with </>). It is perfectly 
clear that, during the evolution ending with the stationary state, the sandpile grows 
under the action of the source, so that it is not reasonable to require that u = <f) in 
the part of the boundary not reached by those transport rays along which no sand is 
poured. For this reason we relax the boundary condition, requiring it < on d£l, and 
by selecting the region r f C d£l where u = (f) in terms of / and of the geometry of the 
transport rays. 

The region Tf also dictates the test functions in the weak formulation of the continuity 
equation. Namely, we require that a solution (it, v) to ([I]) has to satisfy 



v(Du, Di>) dx= fipdx, G C~ (R n \ T f ) , 

n Jn 

thus taking into account the fact that the sand cannot exit from d£l \ Tf. 

We present all the results in a more general setting, which takes into account the 
possibility of homogeneous anisotropies (see also [HI El QUI E] ) • More precisely, we shall 
consider the following system of PDEs in the open, bounded and connected set Q C ]R n 
with Lipschitz boundary: 



(2) 



- div(v Dp(Du)) = f in ft, 

p(Du) < 1, v > in ft, 

(l-p(Du))v = inn, 

u < (p on dfl, 

u = 4> on r^, 



in the unknowns v G L 1 (0), u G W 1,0O (Q,). (Here and in the following we understand 
that u G W 1,00 (Q) denotes the Lipschitz extension to Q of u.) In this formulation: 

- p: W 1 — > [0, +oo) is the gauge function of a compact convex set K C W 1 , of class 
C 1 and containing the origin in its interior; 

- / G ^(O), / > 0; 

- (/): d£l —> M + is a lower semicontinuous function, cf) ^ +oo. 

The plan of the paper is the following. After recalling some notation an basic results, 
in Section [3] we give the precise formulation of the problem, showing that under the 
assumptions listed above we do not have, in general, neither existence nor uniqueness 



of solutions (see Examples 3.4 and 3.5). In order to overcome these obstructions we 
then introduce an additional geometric assumption (see (H5) below), which guarantees 
that the mass is transported straight to the boundary. This condition is automatically 
satisfied if <fi = 0, while in the general case may fail, possibly causing both concentration 
of mass transportation on sets of lower dimension (described by measure-type transport 
densities) or branching of transport paths (and thus multiplicity of transport densities). 
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In Section [4] we prove that ^ always admit solutions, showing that there exists a non- 
negative function vj G such that the pair (u^,,Vf) is a solution (being the Lax- 
Hopf function defined in ^ below) . The main ingredient for this step is a disintegration 
formula for the Lebesgue measure proved by S. Bianchini in [3j. In addition, we prove 
a preliminary but fundamental uniqueness result: if (it<£,i?) is a solution to ([2]), then 
v = Vf. In order to get this result, we show that we need to strengthen the geometric 
assumption (H5) (see Example 3.8 and assumption (H6) below). 

Section [5] is devoted to the characterization of all solutions and, consequently, to the 
uniqueness result. More precisely, we show that there exists a minimal profile Uf such 
that every solution to ^ is of the form (u,Vf) with Uf < u < u$. In particular, the 
u-component is unique, whereas the u-component is unique (and coincides with u^) if 
and only if the support of the source / covers the set of the endpoints of the transport 
rays. Finally, in Section [6] we briefly rewrite our results in the isotropic case Q, and we 
discuss their interpretation in terms of granular matter models. 



2. Notation and preliminaries 

General notation. The standard scalar product of x, y G M n will be denoted by (x, y), 
while \x\ will denote the Euclidean norm of x. Concerning the segment joining x with 
y, we set 

{x, yj := {tx + (1 - t)y; t G [0, 1]}, jx, y\ := {x, yj \ {x, y}. 

Given a set A C M n , its interior, its closure and its boundary will be denoted by int A, 
A and dA respectively. 

We shall denote by C n and % k respectively the n-dimensional Lebesgue measure and 
the /c-dimensional Hausdorff measure. Given a measure p and a //-measurable set F, the 
symbol p[F will denote the restriction of p to the set F. 

If (7 : O — > K is a measurable function, we shall denote by spt g the essential support 
of g, that is the complement in fi of the union of all relatively open subsets 4 C such 
that g = a.e. in A. Notice that spt g is a relatively closed set in f2, but need not to be 
closed as a subset of W 1 . 

Convex geometry. Let us now fix the notation and the basic results concerning the 
convex set which plays the role of gradient constraint for the u-component in ^ . In the 
following we shall assume that 

(3) K is a compact, convex subset of M. n of class C , with G int K. 

Let us denote by K° the polar set of K, that is 

K° := {p G R n ; (p, x) < 1 Vx G K) . 

We recall that, if K satisfies then K° is a compact, strictly convex subset of W 1 
containing the origin in its interior, and K 00 = (if ) = K (see, e.g., |19j). 
The gauge function p : M. n — > K of K is defined by 

p{£) := inf{t > 0; f G tK} = max{(£, n), n G K } , £ G M n . 

It is straightforward to see that p is a positively 1-homogeneous convex function such 
that K = {£ G E" ' : p{C) < 1}. The gauge function of the set K° will be denoted by p°. 



A BVP IN GRANULAR MATTER THEORY 



5 



The properties of the gauge functions needed in the paper are collected in the following 
theorem. 

Theorem 2.1. Assume that K C W 1 satisfies |5|). Then the following hold: 

(i) p is continuously differ entiable in W 1 \ {0}, and 

p(z + v)<p(0 + p(v) ve,r, £ r. 

(ii) K° is strictly convex, and 

P (d + V)<p (0+P°(ri), V^eT 
P°(Z + V) = P°iO + P°{v) 3A>0 : £ = Ar/ or rj = A£. 
(Hi) For every £ ^ 0, -D/?(£) belongs to dK°, and 

(Dp(0, = P (0, (p, < P(0 Vp £ if , p + Dp(0 . 
Proof. See [E], Section 1.7. □ 

In what follows we shall consider M n endowed with the possibly asymmetric norm 

'ii), the unit ball K° of p° is strictly convex but, 
under the sole assumption (|3j), it need not be differentiable. Moreover, the Minkowski 
structure (M n , p°) is not a metric space in the usual sense, since p° need not be symmetric 
(for an introduction to non-symmetric metrics see [II])- Finally, since K° is compact 
and £ int/f°, then the convex metric is equivalent to the Euclidean one, that is there 
exist ci, C2 > such that ci|£| < p°(£,) < C2|£| for every £ £ M n . 

Curves. In the following will denote an open, bounded, connected subset of W 1 with 
Lipschitz boundary. Let us denote by T V)X the family of absolutely continuous paths in 
0, connecting y to x: 

T y>x := {7 G AC([0, i],H), 7 (0) = y, 7 (1) = x} . 

For every absolutely continuous curve 7: [0,1] — > M. n , let us denote by £(7) its length 
with respect to the convex metric associated to p°, that is 

L( 7 ) := f p\i(t))dt. 
Jo 

Since £1 is a compact subset of M n , by a standard compactness argument we have 
that for every x, y £ £1 there exists a (distance) minimizing curve 7 £ Ty,x such that 
£(7) < £(7) for every 7 G r y>a . (see e.g. HJ Thm. 4.3.2], [7, §14.1]). 

The main motivation for introducing the convex metric associated to p° is the fact 
that the Sobolev functions with the gradient constrained to belong to K are the locally 1- 
Lipschitz functions with respect to /?°, as stated in the following result (see [TBI Chap. 5]). 

Lemma 2.2. Assume that the set K C R n satisfies ^ . Let p° be the gauge function of 
K° , let ft C R n be a Lipschitz domain, and let u: Q — > K. Then the following properties 
are equivalent. 

(i) u is a locally 1- Lipschitz function with respect to p°, i.e. 

(4) ^(^2) — u(x\) < p°(x2 — x\) for every \xi,xz\ C fl 

(ii) u £ W 1 ' 00 ^), and Du(x) £ K for a.e. x £ £1. 

(hi) u(x) — u(y) < £(7) for every x,y £ $7 and every 7 £ T y ^ x . 



p°(x — y), x, y £ W 1 . By Theorem 



2.1 
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3. Formulation of the problem 
In this section we shall give the definition of solution (u, v) to the PDEs system 

- div(v Dp(Du)) = f in ft, 

p(Du) < 1, v > in O, 

(5) { (1 - p(Du))v = in ft, 

u < (ft on 3ft, 

u = (ft onT f. 

The basic assumptions are: 

(HI) ft is an open, bounded, connected subset of M. n with Lipschitz boundary; 
(H2) p is the gauge function of a convex set K C W 1 satisfying Q; 
(H3) / belongs to L;L(ft), the set of non-negative integrable functions in ft; 
(H4) (ft: dft — > (-co, +oo] is a lower semicontinuous (l.s.c.) function, (ft =k +oo. 



As we shall see in Examples 3.4, 3.5 and 3.8, this set of assumptions is not enough in 



order to have existence and uniqueness of solutions. Two additional assumptions on the 
geometry of the problem will be introduced in the remaining part of this section. 

The set r j C 9ft where the u component is forced to agree with the boundary datum (ft, 
is dictated by the data of the problem. Using the terminology of the Optimal Transport 
Theory, T f corresponds to the set of initial points of those transport rays on which the 
transport is active. For a rigorous definition, some additional notation is in order. 

Let : ft — > M be the Lax-Hopf function defined by 

(6) u^(x) := inf {(ft(y) + L( 7 ) : y £ Oft, 7 £ T^} , x £ H. 

It is clear that for every the infimum in ^ is attained, that is there exist y G Oft 

and a minimizing curve 7 G ^y,x such that u^(x) = (ft(y) + -^(7). 

Definition 3.1. For x G f2 we call geodesic through x any curve 7 G ^y,x, V £ dQ, 
satisfying u^ix) = (ft(y) + L(j). Moreover, we say that a geodesic through x is (forward) 
maximal if its image is not a proper subset of the image of another geodesic through x. 

We recall that u$ is a Lipschitz function in S7, p(Du ( f l ) = 1 a.e. in VL, and it is the 
maximal function in the space X defined by 

X := {u£ W 1 ' 00 ^) : Du G K a.e. in O, u < (ft on dtl} . 

(Since Q, has a Lipschitz boundary, we understand that all functions in W 1,00 (Q) are 
extended to Lipschitz continuous functions in Q.) We shall show that the Lax-Hopf 
function is always an admissible u-component in ([5]), regardless of the source / (see 



Theorem 4.1 ) 



The function dictates the geometry of the transportation, that is the mass produced 
by the source / runs along the geodesies associated to u^, and falls down at their initial 
points. 

For every ifUwe denote by n(x) the full set of projections of x on dVt, i.e. 
(7) n(x) := {y G 30 : 3 a geodesic 7 G T y ^ x through x} , 
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whereas H M (x) will denote the set of maximal projections, i.e. 

II M (x) := {y G dil : 3 a maximal geodesic 7 G through x} . 
Let us define the set of the initial points of maximal geodesies 

(8) r<£ := {y G dfl: 3 a maximal geodesic 7 G x G $7} = [J TL M (x) 

and the subset of T^ of the initial points of those maximal geodesies where the source / 
is active 

(9) r j := {y G <9f2: 3 maximal geodesic 7 G T y>Xl x G spt /} = n M (x) . 

The set T f turns out to be, in the sandpile problem, the actual portion of (depending 
on the source of matter /) where the sand falls down. The set dQ, \ is the part of the 
boundary closed by the walls. 

We are now in a position to fix the rigorous meaning of problem ([5]). The functional 
setting for the unknowns (u,v) in ^ is Xf x L+(Q), where 

X f := {u G X : u = (j) on T f } , L\(Q) := {v G L 1 ^) : v > a.e. in U} . 

Definition 3.2. A pair (n, v) is a solution to §5§ if 

(i) (u,v)GX f xL\(n); 

(ii) (1 — p{Du))v = a.e. in fi; 
(Hi) for every V G (R n \ F/) 



v(Dp(Du), Dip)dx = f 
Jn 



ftp dx 



Remark 3.3. This notion of solution generalizes the one given in [13] for the table problem 
with walls. The weak formulation in (iii) corresponds to the continuity equation for the 
optimal transport problem subjected to the condition that the mass can flow away from 
Q only through Tf. 

The following two examples illustrate some points that should be taken into account 
in order to deal with existence and uniqueness results for ([5]). In order not to interrupt 
the main flow of the exposition, the details are postponed to Section [7j 

Example 3.4. Under our general assumptions (H1)-(H4), problem ([5]) need not admit a 
solution in the sense of Definition |3.2| Namely, let f2 := fii U ^ Cl 2 , where 



Sll := |(rcos6», rsintf) :0<6><^, 0<r<l}, 
Q 2 ■= [ x = (xi,x 2 ) el 2 : -1 < X\ < 0, \x 2 \ < 1} 
(see Figure [I]) . Let (j) : dQ — > M. be the function 



<t>{x) 



0, x G S := [-1,0] x {-1}, 
+00, otherwise. 



s 
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Figure 1 . The set O in Example [3^4 



Let us consider the isotropic case (p(£) 
Lax-Hopf function is 




|£|) with a constant source / = 1. The 



if x G Oi, 
if x G 172, 



so that = T j = (—1,0] x{— 1}, and the geometry of the geodesies is the one depicted 
in Figure [T] right. In particular the mass collected in the region Oi is transported to 
the origin, and from the origin there is a unique transport ray R going to T^. The 
concentration of the mass on the single ray R corresponds to the fact that the transport 
density v should be a measure with a singular part concentrated on R, which is not 
allowed in Definition 13,21 



Example 3.5. Even in the case of existence of solutions to ([5]), the uniqueness of the 
transport density v may fail if geodesies can bifurcate in the interior. Namely, let 0, := 
Ui U n 2 U f2 3 C M 2 , where 

p2 



(see Figure [2]) . Let < 



:= {x = (xi,x 2 ) G 

-{ 
-{ 

>: dn 



7T 



rcost/, rsm 



rcostf, rsm 



-1 < x\ < 0, 
<9<l, 



x 2 \ < 1} , 
< r < l} , 

be the function defined by 



7T 



< e < 



d>(x) 



0, x G Si U S 3 , 
+oo, otherwise, 



where Si := [-1,0] x {-1} and S 3 := {(cos 6, sin 9) : -vr/2 < 9 < -vr/4}. Let us 
consider the isotropic case (p(£) = |£|) with a constant source / = 1. 

In this example the multiplicity of -u-components depends on the fact that geodesies 
are not forward unique. Namely, given a point y G S3 and a point z £ Q 2 , the curve 
7 := ly, 0] U [0, z\ is a geodesic. It is then clear that we have a lot of geodesies branching 
at 0, so that the mass collected in the region Q 2 is transported to the origin, and from 
the origin it can be distributed in infinitely many ways to any point of S3. 
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ft 2 

► 




^3 

S 3 


Si 






Figure 2. The set ft in Example 3.5 



In order to exclude the phenomena depicted in Examples 3.4 and 3.5, we propose the 
following additional geometric assumption: 

(H5) For every x G ft and for every y G Ii{x) the segment ]y, x\ is contained in ft. 

As a consequence of (H5), every geodesic through a point x G ft is maximal and 
its support is a segment. In particular, the geodesies cannot bifurcate (i.e., they are 
non-branching in the interior). 

Under the assumption (H5), the Lax-Hopf function can be written as 

(10) u^{x) := min {(p(y) + p°(x - y) : y £ dCl] , x G fl, 

while 

= {y G dft : 3 x G ft s.t. fy, x{ C ft and u^(x) = 4>(y) + p {x — y)} , 
Tf = {y£ dft : 3 x G spt / s.t. ]y, i[c(! and u ( f ) (x) = <p(y) + p°(x — y)} . 

For every x G ft, let A(x) be the set of directions through x 

x-y 



A(x) :- 



p°(x - y) 



: y G U(x) 



x G ft . 



where H(x) is the set of projections defined in ([T]) that, under assumption (H5), can be 
written as 

U(x) = {y G T : U(j> {x) = <j>{y) + p\x - y)} . 

Let D C ft be the set of those points with multiple projections, that is 

D := {x G ft : A (a;) is not a singleton}, 

and for every x G ft\D, let p(x) and d(x) denote the unique elements in U(x) and A(x) 
respectively, i.e. 



(11) 



{p(x)} = U(x), {d(x)} = A(x), xeft\D. 
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It can be easily checked that grows linearly along every segment joining x E to 
y E H(x). Let us denote by b(x) be the normal distance from the set D, defined by 



b{x) 



J sup{i > 0; u ( f ) (x + sd(x))) = u ( p(x) + s, Vs E [0, t]}, x E £1 \ D , 

jo xeD, 

and let J be the set 

J ■= U q(%), q(x) ■= x + b(x)d(x), 
where we understand that q(x) = x if x E D. 

Definition 3.6 (Transport ray). We shall call transport ray through x E O any seg- 
ment \p, q(x)J, p E U(x). If \p,q} is a transport ray, the points p and q will be called 
respectively the initial and the fmai point of the ray. 

It is clear from the definition that, if x E O \ D, then there is a unique transport ray 
[p(x),g(x)J through x. In this case there exists a unique number a{x) E (— oo,0) such 
that 

p(x) = x + a(x) d(x) {x E Q \ D). 



Moreover, if assumption (H5) holds, according to Definition 3.1 the segment |p(x),x] is 
a maximal geodesic through x. On the other hand, if x E D, any segment [p, xj, with 
p E H(x), is a transport ray through x. 

The transport rays correspond to the segments where grows linearly with maximal 
slope, i.e. 

(12) u^x + t d(x)) = u^x) + t \/x E Q\D, Vi E [a(x),b(x)] . 

Let £ be the set of those points where is not differentiable. The relationships 
between the singular sets related to the problem are the following (see [3] and |12| 
Prop. 6.4]). 

Proposition 3.7. The sets S, D and J have zero Lebesgue measure. Moreover, DcS 
and D C J C D, with possibly strict inclusions. In addition, if ua is differentiable at 
if!!, then x has a unique projection and A(x) = {Dp(Du ( j ) (x))} . 

Under the assumptions (H1)-(H5) we will be able to construct a mass transport 
density Vf such that the pair (u$,Vf) solves the system ^ (see Section[4|. Unfortunately, 
these assumptions are not enough to ensure the uniqueness of the mass density v, due 
to the possibility of transporting a fictitious amount of mass along those transport rays 
}p,q{ with bot h endpoints on d£l. The weak formulation of the continuity equation 



(Definition 3.2 iii) ) prevents this possibility only for q E dfl\Tf. The following example 
shows that a single q E Tf may be the final point of a set of rays covering a region with 
positive Lebesgue measure. 

Example 3.8. Let O := fii U VL 2 U fi 3 C M 2 , where 

«1 := (-1,0) x (-1,1), ft 2 := [0,1) x (0,1), 
n 3 := I (r cos 6,r sin 6) : r E (0,1), ~ < 6 < -^j , 
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and let eft: d£l ->Mbea function such that, in the isotropic case with source / = 1, the 
geometry of transport rays is the one depicted in Figure [3j In particular = Tf = 
Si U S2 , where 

Si := ((-1,0] x {-l,l})U{(cos0,sin0) : -| < < -J}. 
S 2 :=((0,l)x{0})U({l}x(0,l)) 



(see the details in Section[7j). Notice that (0, 0) S is the final point of all the transport 
rays covering O3. Hence, we can construct a w G L^Os) such that div(w; Du^) = 0, 
which can be added, once prolonged to zero in f2 \ O3, to any admissible v, loosing the 
uniqueness of mass transport density. 



In order to exclude the behaviour described in Example 3.8 we need the following 
additional assumption: 

(H6) (Jndti) n I> = 0. 
It is worth to remark that, as a consequence of (H5), one has (J n d£l) n = (see 
Lemma 3.9 below), so that (H6) can be viewed as a mild additional assumption in order 
to separate the sets J n dfl and I"^. 

Lemma 3.9. Assume that (HI), (H2), (H4) and (H5) hold, and let {pi,qi}, \p2,Q2\ be 
two distinct non-trivial transport rays. Then q\ ^ P2- In other words, J D = 0. 



Proof. Assume by contradiction that q\ = 
support is [pi,p 2 ] U [p2,92]- By (III) we get 



P2, and let 7 6 r pii92 



be the curve whose 



= U4,(P2) + P u (q2 -P2) = Ucf>(pi) + P V (P2 -Pi) + P u fe -P2) = u^ipi) + L(7), 



A<?2 ~ P2) = Ucf>(Pl) + P°(P2 ~ Pi) + p° 
hence the curve 7 is a geodesic, in contradiction with (H5) 

4. Existence of solutions 



□ 



For the reader's convenience we collect here all the assumptions that have been intro- 
duced in the previous section. 
(HI) f2 is an open, bounded, connected subset of W 1 with Lipschitz boundary; 
(H2) p is the gauge function of a convex set K C R n satisfying 
(H3) / belongs to L\(Q)) 

(H4) (j>: dCl — > (—00, +00] is a l.s.c. function, <f> ^ +00; 

(H5) For every x £ Q and for every y S U(x) the segment \y, x\ is contained in f2; 
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(H6) (Jnffi)nr^ = 0. 

This section will be devoted to the proof of the following result. 

Theorem 4.1. Let (Hi)-(H6) hold. Then there exists a unique Vf G Li_(0) such that 
the pair (xtfaVf) is a solution to the PDEs system 

This theorem gives a partial uniqueness result for the transport density, since it states 
that the u-component associated to the profile u = is unique. We shall see in SectionlB] 



Theorem 5.5 that, indeed, Vf is the only admissible w-component of every solution (u, v) 



Since G Xj, and p(Du^) = 1 a.e. in Q, in order to prove the existence part of 



Theorem 4.1 it is enough to show that 



(13) / v } (Dp(D U(j> ), Dili) dx= fipdx, G C c °° (R n \T f ). 

Jn Jn 

This can be done with a minor effort, since it turns out that we can choose as v- 

component the same function constructed in [3l Q2] for the case of admissible boundary 

data, that is the one that satisfies 



(14) / v f (Dp(Du$), D^)dx = / fipdx, V^eC c °°(ft). 

Jn Jn 

In order to write the explicit form of the function Vf, we need the fundamental result 

concerning the disintegration of the Lebesgue measure along the transport rays due to 

S. Bianchini (see (3j Thm. 5.8]). The proof of this formula is mainly based on the fact 



that the divergence of the vector field d defined in ( 1 1 ) is a locally finite Radon measure 
in VL. Moreover, decomposing divd into its absolutely continuous and singular part 
(w.r.t. the Lebesgue measure), divd = (divd) ac £ n + (divd) s , it turns out that (divd) s 
is a positive measure in f2 (see [3], Section 5). These properties allow to describe the 
evolution of the Hausdorff measure of the (n — l)-dimensional sections of "cylinders" of 



transport rays (see Remark 4.3). 



Theorem 4.2. Let (HI), (H2) and (HA) hold. Then there exists a sequence {A^}^^ 
of subsets of Q with the following properties. 

(i) A k C (fi \ D) n {{x, ej k ) = Zfc} for some jk G {1, . . . , n} and G M. 

(ii) A k is measurable w.r.t. the (n — 1)- dimensional Hausdorff measure. 

(iii) The sets T k := \J x€ ^ k }p(x), q(x){, k G N, are pairwise disjoint, Lebesgue mea- 
surable subsets of Q, and C n (Q, \ [j k T/-) = 0. 



(iv) For every h G L 1 ^), the function t h-> h(x + td(x)) a(tdj k (x),x), where d 0k 



.x 



(d(x), ej k ), belongs to L 1 (a(x) , b(x)) for every x G \] k Ak- Moreover, the disin- 
tegration formula 

hdx = Y] I I I ( } h(x + td(x))a(td jk (x),x)dt) d jk (x) dH n ~ 1 (x) 



(15) 



H k J A k \ Ja(x) 

holds, and for T-L n ^ l -a.e. x G U fc A k the function a(-,x) is the solution of the 
linear ODE 
( (i 

—a(td jk (x),x) = (divd) ac (x + td(x))a(td jk (x),x), t G (a(x),b(x)), 
a(0, x) = 1 . 
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Moreover, this solution is strictly positive. 



We remark that the decomposition introduced in Theorem 4.2 is clearly not unique. 
In the following we shall always assume that such a decomposition has been fixed. Once 
the decomposition is given, we shall use the notation 

a(x + td(x)) = a(tdj k (x), x) dj k (x), k G N, x G A k , t G (a(x), b(x)), 

so that a is defined and strictly positive on the set 



(16) 



J)' := I x + id(x) : xe\jA k ; t€ (a(x),b(x)) \ ={jT k 
{ k J k 



i.e., by Theorem |4.2[ iii) , almost everywhere on Q. For every x G f2' the ODE (15) along 
transport rays can be written as 

^ a ( x + td(x)) = (div d) ac (x + td(x)) a(x + td(x)), 

and a = 1 on |J fc ^4fc, while the disintegration formula for the Lebesgue measure given 
in Theorem |4.2[iv) simply becomes 



(17) f hdx = Y] [ (/ ] h(x + td(x))a(x + td(x)) dt] dH^ix). 

JO k ^ A k \Ja(x) J 

Remark 4.3. It can be of interest to understand the geometrical meaning of the function 
a(t,x) in Theorem |4.2| (iv) . For fixed k £N and t € M, let us define the map 

tpt-.Ak^W 1 , x (->• <p\x) := x + t^jK , 

d jk (x) 

and let A\ := cp (A k ). Since A k is contained in the hyperplane {x : (x, e Jfe ) = z k }, we 
have that A\ C {x : (x, ej k ) = z k + i}. Let us consider the measure fi on A k such that 
the push-forward (</?*)jjM of fi is H n [A k , that is, 

(18) fi (((p t y 1 (F)) = % n - l {F) for every U n ~ ^measurable set F C A\. 



It turns out that the measure fi defined by (18) is absolutely continuous w.r.t. 7-L n ~ [A k 
and fi = a{t,x)'K n ^ l \ L A k (see [3j Lemma 5.4]). 
Let C C f2 be a "cylinder" of rays of the form 

C = {x + sg?(x) : x G j4fc, s G (a(x), &(x))} , 

and let xc denote its characteristic function. Then, for every h G L 1 (C), using Fubini's 
theorem and the definition of fi we have that 

hd£ n = [ [ h XC dH n - l dt 



Jo 



(hxc)(<p\x))a(t, x) dH^Hx) dt 

b(x) 

h(x + s d(x))a(sdj k (x), x)dj k (x) ds dl-L n (x) 

a(x) 
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Hence, roughly speaking, the disintegration formula (17) corresponds to covering with 
cylinders of transport rays and applying Fubini's theorem on each cylinder. 

We are now in a position to show that the problem of finding the solution v to the 
equation 

v(Dp(Du (t> ), Dip) dx= I fipdx, Vtp G C c °° (M n \ T f ) 



(19) 

can be solved by the Method of Characteristics. We start recalling another result proved 
in [31 Sect. 7]. 

Theorem 4.4. Let (HI) -(HA) hold. If v G L\(fL) satisfies 



(20) 



v(Dp(Du ( f ) ), Dip) dx 



fipdx, Vip£C™(Q), 



then for a.e. x £ £1 the function v is locally absolutely continuous along the ray t \— > 
x + td(x), t G (a(x), b(x)), and satisfies 

d 



(21) —[v(x + td(x))a(x + td(x)) 



-f(x + td(x)) a(x + td(x)) , t G (a(x), b(x)). 
Moreover, if the final point q(x) := x + b(x)d(x) belongs to ft, then 



(22) 



lim v (x + td(x)) a(x + td(x)) = 0. 
t-^(z)- 



Finally, the function vj: 0, — > R defined by 
(23) v f (x) :-- 



b{x) ,/ u ^ a(x + td(x)) , 

f{x + td(x)) — . . v dt, 

a(x) 



a.e. x G . 



belongs to L^(Q) and is a solution to (20). 



Remark 4.5. The fact that vt satisfies (21) along almost every ray can be easily proved 
observing that, for every x E Q', 



Vf(x + td(x)) 



b(x) 



f(x + sd(x)) 



a(x + sd(x)) 
a(x + td(x)) 



ds, 



so that 



(vf a)(x + td(x)) 



b(x) 



(f a)(x + sd(x)) ds. 



From this last equality we also see that Vf satisfies the terminal condition (22) for every 
x G 17', and not only for those points x £ Q satisfying q(x) G O. 



Remark 4.6. From the definition (23) of Vf and the fact that a(-,x) is strictly positive 
in (a(x),b(x)), we deduce that the essential support of vj coincides with the closure of 
the set 

{U]p(x),x[: xGspt/}. 
In particular Vf = along all the transport rays starting from the points of \ Tf. 
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Remark 4.7. Notice that neither (H5) nor (H6) are needed in order to find a solution to 



(20). Namely, the weak formulation (|20|) corresponds to the continuity equation of the 



optimal mass transport problem in which the mass can freely flow away from Q, at the 



first time it touches the boundary. Hence, the phenomena depicted in Examples 3.4 and 



3.5 cannot happen. 



Since all functions tp £ C£°(fi) (extended to in M. n \ Q) are admissible in the weak 
formulation ( 19 ) of the transport equation, it is clear that every solution v to ( 19 ) satisfies 



(21) along almost every ray. 

Using the disintegration formula for the Lebesgue measure it is not difficult to prove 



that the function vj defined in (23) above is a solution to (13), provided that the mass 
can flow away at the initial points of the transport rays, i.e. 

L(x): xe\jA k , ]p(x),g(*)[nspt(/)^0j QT f . 

This condition is clearly satisfied if (H5) holds. 



Theorem 4.8. Assume that (H1)-(H5) hold. Then (u<f,,Vf) is a solution to (14). 



Proof. From Theorem 4.4, the function vj satisfies the ODE (21 ) along almost every ray 



x + td(x), t 6 (a(x), b(x)). Moreover, as already observed in Remark |4.5 



(24) (v f a)(q(x)) := lim (v f a)(x + td(x)) = 0, 

t— >b(x)~ 



%" -1 -a.e. x 



e\jA k . 



n' defined in ( 16 ). 



It is not restrictive to assume that (21) and (24) hold for every x belonging to the set 



3.7 



£1' C £1 \ J C 0, \ D, we have that, for every x G 0,' , 



Since, by Proposition 
d(x) = Dp{D U4> {x)) and, for" V e Cf (R n \T f ), 

{Dp{Du 4> )){x + td(x)), Dip(x + td(x))) = (d(x), Dtp{x + td(x))) = rp'{x + td(x)) , 



where the prime denotes differentiation w.r.t. t. Hence the disintegration formula (17) 
gives 

/ Vf(Dp(Du ( f ) ), Di/j) dx 
Jn 

( 25 ) r rb(x) 

= J2 / {v f a){x + td{x))^{x + td{x))dtdrl n - l {x). 
k J A k Ja(x) 



An integration by parts in the inner integral of (25) and formula (21) lead to 

/ (vfa)(x + td(x))ip'(x + td(x)) dt 
( 26 ) - H x) 



(vfaif))(q(x)) — (vfai[))(p(x)) + / (fa)(x + td(x))ip(x + td(x)) dt , 

J a(x) 



where we use the convention 



(vfa)(p(x)) := lim (v ja)(x + td(x)). 

t— >a(x)+ 
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If p{x) belongs to Tf then the test function ip vanishes in a neighbourhood of p(x). On 



the other hand, if p{x) Tf, then / vanishes on the ray g(x)[, so that, by (21), 

Vfa is constant along that ray. Since (v fa)(q(x)) = 0, then also {vfa)(p{x)) = 0. 

Finally, using again the disintegration formula, the weak formulation (19), (24), and 
the boundary conditions (vfaip)(p(x)) = (vfa)(q(x)) = 0, we conclude that Vf is a 
solution to (19), i.e. the pair (u,p,Vf) is a solution to §5§. □ 



In order to prove a uniqueness result we need to exploit assumption (H6). 



Theorem 4.9. Assume that (H1)-(H6) hold. Then a function v G L\(ft) satisfies (19) 



if and only if, for a.e. x £ Vt \ D , the function v is locally absolutely continuous along 
the ray 1 1-)- x + td(x), t G (a(x), b(x)) and 

< d 



(27) 



[v(x + td(x)) a(x + td(x))] = - 

lim v(x + td(x)) a(x + td{x)) 
t-±b{x)- 



-f(x + td{x)) a(x + td(x)), 
= 0. 



As a consequence, the function vj defined in (23) is the unique solution of (19). 



Proof. By Theorem [Ol hi) it is clear that two functions both satisfying (27) along almost 



every ray must coincide almost everywhere. Hence the uniqueness result will be achieved 



once we prove that every solution to (19) satisfies (27). 
Let v € L\(p) be any solution to (19), and let ip £ C c c 



\I7). Notice that v is a 



solution to ( |20| ); by Theorem 4.4, it is not restrictive to assume that ( |21[ ) holds for every 
x belonging to the set defined in (16), whereas (22) holds for every x G Q' such that 
q(x) G fl 



Reasoning as in the proof of Theorem 4.8 (see (25) and (26) with v instead of Vf), 
and using the fact that v is a solution to (19]), we obtain 



E 

k 



[(vaip)(q(x)) - (van/))(p(x))] dn n -\x 



0. 



By (22), if q{x) belongs to O, then (va)(q(x)) = 0. Moreover, if p(x) belongs to Tf then 
the test function -0 vanishes in a neighbourhood of p(x). Hence the formula above is 
equivalent to 



(28) 



E 

k 



(va^)(q(x))dH n ' l (x) 



E 

k 



{va^){p{x))d% n - l {x) = 



where 

B k := {x G Ak : q(x) G dtt} , C k := {x G A k : p(x) Tf} 
Let us consider the following subsets of dfl: 



B 



L(x): xGljBfej C JndQ, C:=L(x): x G J C k \ C T \ T j 



From the very definition of C, (H6) and the inclusion Tj CT^we deduce that 

(29) Fnr / = 0, Hnc = 0. 
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Hence, choosing 5 > such that B$(B) n C = B$(B) n T j = 0, by (29) we can construct 
a function tp G C£°(M n ) satisfying 

^>0 onB, tp = on \R n \B s (B)]DT f \JC, 

so that tp G (M n Recalling that ua > 0, and using -0 as a test function in (28) 

we conclude that (va)(q(x)) = for % n ~ -a.e. x G G N, so that the boundary 

condition in d2Tb holds. 



□ 

Another easy consequence of both the representation formula for Vf and the disinte- 
gration (17) is the following stability result. 

Theorem 4.10 (Stability). For every /i,/2 G Li(fi) i/iere Zio/rfs 
IKi ~ VfallLHtt) < diam(O) ||/i - / 2 || L i (f7) , 
where diam(fi) denotes the diameter of the set Q. 



Proof. From the disintegration formula (17), the representation formula (23) and Re- 
mark |4?5] we have that 

rb(x) \ 

\(v h - v h )(x + td(x))\ a(x + td(x)) dt ) dU n ~ l {x) 

b(x) \ 




diam(f2) 



(|/i — /2I a)(x + sd(x))ds J dt 

b(x) 



dn n -Hx) 

(l/i - h\ a)(x + sd(x))ds ] c#T _1 (x) 



= diam(fi) ||/i - / 2 ||xi(n) > 
completing the proof. 

5. Uniqueness of the solutions 



□ 



Up to now we have proved that there exists a unique v = v f such that the pair {u^,vA 
solves ([5]). In this section we shall prove that, actually, Vf is the unique admissible v- 
component for ([5]). 

For what concerns the u component, we start by proving that the Monge-Kantorovich 
system §5§ is the Euler-Lagrange condition for the minimum problem with gradient 
constraint 



(30) 



mm 



/— J fudx: u G Xf\ 



It is clear that, since / is non-negative and is the maximal element in Xf, then 
is a solution to (30). Moreover every u which minimizes (30) has to agree with on 
spt/. 



Theorem 5.1. Assume that (H1)-(H5) hold. The minimum problem (30) and the 
system of PDEs §5§ are equivalent in the following sense. 

(i) u G Xj is a solution to (30) if and only if there exists v G L+(£l) such that (u,v) 
is a solution to 
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(ii) Let u G Xf be a solution to (30). Then (u,v) is a solution to ^ if and only if 
(u^jv) is a solution to Q. 

Proof. The proof is similar to the one of Theorem 5.3 in [12]; for the reader's convenience 
we sketch here the main steps. 

In what follows we will freely use that fact that, by a density argument, the difference 
u — w of two functions u, w G Xf can be used as test function in the weak formulation 
(19) of the transport equation. 

Let us denote by Ik the indicator function of the set K, that is 

'o iixeK, 



I K {x) ■ 

so that 



+oo if x G W 1 \ K . 



F(u) := [ [I K {Du) -fu)dx = - [ fudx, Vu G X f . 
Jn Jn 
Since the gauge function p is differentiable in IR n \ {0}, the subgradient of Ik can be 
explicitly computed, obtaining 

'{aDp(£):a>0} if £ G dK, 
{0} if^Gintif 

(see e.g. PS Sect. 23]). 

Hence, if (u, v) is a solution to ([5]), by condition (ii) in Definition 3.2 we have that 
v(x)D p(Du{x)) G dlK{Du(x)) for a.e. x G f2, so that, for every w G Xf 



F(w) - F(u) > / v{Dp(Du), Dw - Du) dx - 

Jn Jn 



f (w — u) dx = , 



v f (D p(Du ( p) , Du — Du<f,} dx = — Vf(l — (Dp(Du ( j ) ), Du)) 

Jn 



where the last equality follows from the fact that w — u can be taken as test function 
in (19). This proves that u is a solution to (30). 

Assume now that u G Xf is a minimizer for F, so that f{u — u^) = a.e. in 0, due 
to the maximality of u^ in Xf and the fact that / > 0. Again we can choose u — u^ as 
test function in the transport equation ( [l9| ) solved by u^,, getting 

dx . 

On the other hand, by Theorem |2.1[ ii) and the fact that Du G K a.e. in Q, we have 

1 - {Dp(Du (f> (x)), Du{x)) > 0, a.e. iGfi, 

1 - (Dp(D U(/ ,(x)), Du{x)) = Dp(D U(j> (x)) = Dp(Du(x)), 

so that VfDp(Du ( p) = VfDp(Du) and Vf(l — p{Du)) = a.e. in O, that is {u,Vf) is a 
solution of ([5]) . This concludes the proof of (i) . 

Let us prove (ii). The previous computation shows that if {u^^v) is a solution of 
and u G Xf is a solution of the minimum problem (30), then also (u, v) is a solution of 
([5]). Finally, let (u, v ) be a solution to ([5]). Upon observing that v p(Du) = v a.e. in 
and choosing u — u^ as test function in the weak formulation of the transport equation 

-div(vDp(Du)) = f, inn, 
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we conclude that (u^, v) is a solution to §5§. □ 

As a consequence of the previous results we obtain the following uniqueness result for 
the v component. 

Corollary 5.2. Assume that (H1)-(H6) hold. If(u,v) is a solution to then v = Vf 
in Q, and u = on spt /. 

Proof. Let (u,v) be a solution to ([5]). By Theorem |5.1[ i) we have that u is a solution to 
(30) hence, by Theorem 5.1 ii), (u^,t>) is a solution to ([5]). The conclusion now follows 
from Theorem 14.91 □ 



We now introduce another element of Xf, which will play the role of the minimal 
admissible profile. This minimal profile uj depends on the source / and is defined by 

(31) Uf(x) := sup{u^(z) - L(j) : z G spt(/), 7 G T x>z } , 

with the convention uf = —00 if / = 0. (We recall that, by definition, spt(/) is a 
relatively closed set in Q.) 

Proposition 5.3. Assume that (H1)-(H4) hold. If f then the function Uf defined 
in (31 ) belongs to Xf and Uf = on spt /. Moreover, every function u G Xf such that 
u = Uf in spt(/) satisfies Uf < u < 11$ in tt. 

Proof. By the very definition of Uf, we have that Uf satisfies Q in ft, hence, by 
Lemma |2^2l Uf G W 1 ' 00 ^) and Du f G K a.e. in fl. Moreover, for every 2 G spt(/) 



and y G 9TF we have that 

^>{y)>u <j) (y)>u^(z)-L(i), V7 G Y y>z , 



so that </>(y) > Uf(y). In order to prove that Uf = 4>onTf,lety£Tf and let 2 G spt(/) 
be such that there exists a maximal geodesies 7 G T^, i.e. 1^(2) = (f>(y) + -^(7)- Then 

</>(y) > u f {y) > u^z) - L(i) = 4>{y). 

It remains to prove that, if u G Xf coincides with Uf on spt(/), then u > Uf. (The 
inequality u < is trivially satisfied by the maximality of uj, in Xf.) Namely, for every 

x G ft there exist 2 G spt(/) and 7 G T XjZ such that u/(x) = it^(z) — £(7), hence by 
Lemma |2.2riii) 



n/(x) + L(7) = u${z) = n(2) < u(x) + L(j) 
i.e. Uf(x) < u{x). □ 



Remark 5.4. Since u/ = in spt(/), exploiting the explicit representation formula (23) 
of Vf, we can infer that Uf = in spt(u/) (see also Remark 4.6). 

The following result is the analogous of Theorem 7.2 in [12J. 

Theorem 5.5 (Uniqueness). Assume that (H1)-(H6) hold. Then a function u G Xf is 
a solution to (30) if and only if Uf < u < u^. Moreover, the function Uf coincides with 
in 0, (and hence is the unique solution to (30)) if and only if J C spt(/). 
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Proof. Let us first consider the trivial case / = 0. In this case Xf = X , Uf = — oo and 



every function u £ X is a solution to the minimum problem ( 30 ) . 



Let now assume that / ^ 0. The first assertion follows from Proposition |5.3| and from 



the fact that u S Xf is a solution to (30) if and only if u = on spt(/). 
Let us prove the uniqueness result. 



Let J C spt(/). From Proposition 5.3 we have that uj = on spt(/), hence on J. 
Let x £ VL \ J be given, and let q(x) E J be the endpoint of the ray through x. We have 

u^x) = u <t> {q{x)) - p°(q(x) - x) = u f (q(x)) - p°(q(x) - x) < Uf(x) < u^{x) , 

and hence uj(x) = u ( j ) (x). 

Assume now that = Uf in f2, and assume, by contradiction, that there exists a 
point xq £ J, xq spt(/). By definition, there exist z S spt(/) and a curve 7 £ T XQ)Z 
such that 

(32) u f (x ) = u^(z) - L(7). 



Moreover, by (10), if yo £ n(a?o) one has 



(33) u<t>(xo) = Uj,(y ) + p°(x - y ) 



Since, by assumption, u^Xq) = u/(xq), (32) and p3j) yield 



P°{z - yo) < p°{x - y ) + L(i) = u^z) - u^yo) < p°(z - y ), 
i.e., p°(xq — yo) + -^(7) = P°{z — Uo)- From Theorem |2.l[ii) it follows that [yo, xq} U 7 



[yo, z}, that is ^(7) = p°(z — xq) and xq G [yo, zj. Finally, since z 7^ xq, by (32) and the 



definition of final point we have that xq J, a contradiction. □ 
For the reader's convenience, we summarize here the results we have obtained. 

Corollary 5.6. Assume that (H1)-(H6) hold. Then 

(i) If (u,v) is a solution to (51), then v = Vf. 

(ii) (u, Vf) is a solution to (15) if and only if u £ X and uj < u < in Q. 

(iii) (u<f,,Vf) is the unique solution to ^ if and only if J C spt/. 

6. The isotropic case: application to sandpiles 

In this section we shall focus our attention on the isotropic case p(£) = |£|, £ G M n , 
(n = 2 in the model problem) translating the results of the previous sections in terms of 
description of the equilibrium configurations for sandpiles on a flat table with a vertical 
rim, and comparing these results with the ones known in literature. 

In what follows we always assume that (HI), (H3), (H4), (H5) and (H6) hold. Fur- 
thermore, when considering the sandpile model, one could prefer not to allow profiles 
with negative height. If this is the case, one has to assume that 

(a) The boundary datum (f> (corresponding to the height of the rim) is non-negative; 

(b) The space Xf contains only non-negative functions; 



(c) The minimal function Uf, defined in (31), is replaced by m&x{uf, 0}. 

If the sand is poured by a vertical source / on the table occupying the region Q,, the 
admissible equilibrium configurations can be depicted in terms of the pair (it, v), the 
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profile of the standing layer and the thickness of the rolling layer, which are solutions of 
the PDEs system ([!]), that is 

(i') (u,v)GX f xL\(n); 
(ii') (1 - \Du\)v = a.e. in U; 
(iii') for every ip G C c °° (R n \ T f ) 



v(Du, Dift) dx = / fiftdx, 
n Jo, 

where 

T j = {y £ dO, : 3x G spt / s.t. }y, x{ C $7 and n^(x) = ^(y) + \x — y\} 
is defined in term of the Lax-Hopf function 

u^x) = mm{(ft(y) + \x - y\, y G d£l}. 



Following (10), we have that 

u<j>(x) = min{(/)(y) + \x - y\, y G 1^}, x G Q, 

where 

r (j, = {y G dil : 3x € H s.t. ]y, x[ C and «^,(a;) = 4>{y) + \x — y\}. 

Hence in the sandpiles model, is the set of the initial points of the transport rays where 
the matter is allowed to run down, that is where the u-component could be non-zero. 

On the other hand, is the set of the initial points of the transport rays along which 
the transport is active. It is the effective border of the container, where the standing 
layer fills the gap with the rim of height (ft (since u = Uf = = (ft on T f for every 



admissible profile, see Theorem 5.5) and the exceeding sand falls down. 

Finally, <9f2 \ is the part of the boundary closed by walls that the sand cannot 
overcome, no matter what the source is. 

The results obtained in the previous sections, interpreted in the light of the sandpile 
model, are the following. 

Theorem 6.1. Assume that (HI), (H3)-{H6) hold. Then: 

(i) The thickness of the rolling layer is uniquely determined by the data of the prob- 
lem. 

(ii) The admissible profiles of the standing layer are wedged between the minimal 
profile Uf and the maximal profile u^. In particular, every admissible profile has 
to agree with the maximal one in the region where the transport is active. 

(iii) There exists a unique admissible configuration if and only if the source pours 
sand on the whole ridge of the maximal profile. 

These results extend the known results for the open table problem without walls (see 
[H [5j El [T7]) and with walls (see p3]). Moreover, we have clarified the issue of the lack 
of uniqueness of the rolling layer we have faced in a less general setting in |12j . 
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7. Examples 

In this section we detail some computation that were skipped while presenting Exam- 
ples 3.4| 3.5 and |3.8| in Section [3j 

In all examples the metric is isotropic, i.e. p(£) = |£|. Moreover, we assumed that 
/ = 1, so that the set Tj (defined in ([9])) of initial points of maximal geodesies intersecting 
the support of / coincides with the set (defined in Q) of all initial points of maximal 
geodesies. In the first two examples only assumptions (H1)-(H4) hold, so that maximal 
geodesies need not be segments. In any case, by Proposition 5.3 we have that Uf = u^, 
so that the Lax-Hopf function is the unique candidate for the u-component of the 
solution. 

Concerning the ^-component, by Theorem 4.4 we have that, along almost every trans- 
port ray, the ODE (21 ) is satisfied. Moreover, if the final point q of a transport ray does 
not belong to the closure of the set of initial points, then va must vanish on q (see the 
end of the proof of Theorem 4.9). 

Details on Example 3.4. We have already observed that every solution (u,v) of ^ 
must satisfy u = in f2, whereas v(x) = 1 — X2 for every x = {x\, x-i) E f&2, so that, for 
every test function ip E C^°(M? \ F/), 



(v(D U<j) , Dip) -ip)dx = 0. 



On the other hand, v must be locally absolutely continuous along any ray (0, 1) 3 r i— > 
(r cos 0, r sin 9), < 9 < tt/2, and 



d_ 
dr 



[rv(r cos 9,r sin 9)] 



-r € (0,1), < 9 < 



7T 



(see (|21[)). Then, on Oi, we have that 



1 



uircosf, rsinl 



+ 



c{9) 



re (0,1), < 9 < 



7T 



2r r v ' " 2 

for some integrable function c > 0. Now, an explicit computation in polar coordinates 
gives 

fTr/2 



^{Du^, Dip) — ip) dx 



n 



c(9)ip(cos9,sin9) 



+ c{9) V(0) 



d9 



and this last integral clearly cannot vanish for every choice of the test function ip. 



Details on Example 3.5, The Lax-Hopf function can be easily computed: 



u<t>(x) 



1 + X2, if x E fii, 
1 + |x|, if x E 5^2; 
1 — |x|, if x E O3. 



It is clear that T„ 



(5 1 U5 3 )\{(-l,-l),(2- 1 /2,_ 2 -i/2 )} . 
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v f (x) :-- 



1 + x 2 , if x G fii, 

(l-|x| 2 )/(2|x|), ifxG0 2 , 
Jx|/2, ifxG^ 3 , 



and 
(34) 
with 



w(r cos 6, r sin I 



cGLI (-tt/2, -tt/4), 



c(6)/r on Q3, 
on \ Sl 3 , 



-tt/4 
tt/2 



c(6) d6 



7T 



is a solution to ^ in the sense of Definition 3.2 Notice that, as usual, the function v f 
is obtained solving the ODE (21), while w is a solution to — div(wDu ( f,) = in O, with 
J Q w = £ 2 (0 2 ). Namely, given ip G C£°(IR 2 \T/) and denoting by 



[(vf + w)(Du ( f ) , Dip) — ip] dx, j = 1,2,3, 



passing in polar coordinates we easily get 

i 1 = o, j 2 = -|v(o), i 3 = ^(o) 

hence I\ + I 2 + I3 = and the claim is proved. 



-tt/4 
tt/2 



7T 



c (0)d0 = -V(0) 



Details on Example 3.8, The boundary datum corresponding to the geometry of 
transport rays depicted in Figure [3] is the continuous function <p : dtt — > M satisfying 



4>(x) 





l xeS 2 
1 - y'xf + x\ x £ S3 

X\ 

1 - b2 1 X G S5 



((-1,0] x {-1,1}) U {(cos 9, sin 6) : -f < < -f } 

((0,l)x{0})U({l}x(0,l)) 

{(xi,-xi), xi G (0, ^)} 

(0,1) x{l} 

{-l}x[-l,l]. 



The Lax-Hopf function can be computed as follows: 



it^(x) 



1 — |x 2 | x G Oi 

1 - yjx\ + x\ X G 

x G Q' 2 := {xi G (0, 1), < x 2 < f } 



1 + x 2 

2 -xi x G 2 ' := {x 2 G (0,1), 1 j 

^xf + (i-x 2 ) 2 xen 2 \(n 2 un 2 '). 



(i-M) 



< Xl < 1} 



Hence r^, = = Si U S 2 . Since (H1)-(H5) are satisfied, by Theorem 4.8 the pair 
(u^jVf) is a solution to §5§. 
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On the other hand (0,0) E Tf is the endpoint of all the transport rays covering Q3, 
so that (H6) is not satisfied. Hence, if w : Q — > K is a function of the form (34), with 
c G L\(—ir/2, — vr/4), then for every test function -0 £ C*^ (M ra \ Tf) one has 



w(Du ( j > , Dip)dx = / w(Du,j > , Dip)dx = 0, 



so that the pair (u^^Vf + w) is a. solution to (ph. 
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